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Abstract 

The spin transfer torque (STT) can lead to steady precession of magnetization without any 
external applied field in magnetic spin valve where the magnetic layer have very different spin 
diffusion length. This effect is associated with an unusual angular dependence of the STT, called 
"wavy" (WAD-STT), predicted in the frame of diffusive models of spin transfer. In this article, 
we present a complete experimental characterization of the magnetization dynamics in the pres- 
ence of a WAD-STT. The results are compared to the prediction of the magnetization dynamics 
obtained by single domain magnetic simulations (macrospin approximation). The macrospin simu- 
lations well reproduced the main static and dynamical experimental features (phase diagram, R(I) 
curves, dependence of frequency with current and field) and suggest that the dynamical excita- 
tions observed experimentally are associated with a large angle out-of-plane precession mode. The 
present work validates the diffusive models of the spin transfer and underlines the role of the spin 
accumulation and the spin relaxation effects on the STT. 
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1. INTRODUCTION 



A spin polarized current can exert a torque on the magnetization of a ferromagnetic 
body, via transfer of spin angular momentum. This spin transfer effect, originally introduced 
by J.Slonczewski and L.Berger^l, has opened a new route to manipulate a magnetization 
using an electrical current. Today, it is the subject of extensive experimental and theoretical 
research motivated not only by its fundamental interest but also by promising applications in 
the field of magnetic memories and microwave devices for future telecommunication systems. 
In most experiments, this effect has been studied in pillar shaped F1/NM/F2 trilayers in 
which a magnetic layer Fl with a fixed magnetization is used to prepare a spin polarized 
current injected in a magnetically free magnetic layer F2. When the spin polarized current 
enters the free layer, the part of the spin current that is non-collinear to the magnetization 
of the free layer is absorbed and transferred to the local magnetization; this spin transfer is 
equivalent to a torque exerted on magnetization. In standard structures such as Co/Cu/Co 
or NiFe/Cu/NiFe, at zero or low applied field, the spin transfer torque (STT) leads to an 
irreversible switching of the magnetization from one static configuration to another one^ (so 
called CIMS effect). For larger fields (typically higher than the coercive field of the free layer), 
the spin transfer torque compensates for the damping torque resulting in a steady precession 
of magnetization around the local internal field. Owing to the giant magnetoresistance 
(GMR) effect, this high frequency magnetization precession goes with voltage oscillations 
in the microwave range^. These two behaviors (magnetization switching and precession) 
open two branches of potential technological applications of the spin transfer effect: first, 
a new and more reliable way to write the magnetic bit in magnetic memories and second, 
a new type of sub-micrometric microwave oscillators for future telecommunication devices 
combining high quality factor and large frequency agility with current and field ®. 

Following the original approach of SlonczewskP, several theoretical studies underlined 
that the STT is directly related to the spin polarization of the spin current entering the 
magnetic free layer and more precisely to the absorption of the spin current's component 
transverse to the magnetization^!. In the first generation of spin transfer models^ED, this 
spin polarization was calculated in a ballistic assumption. The magnetic multilayer is con- 
tacted to reservoirs with uniform chemical potentials and the spin polarization of the current 
arises from the spin dependent reflections at the interfaces of the multilayer and from spin 
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dependent scattering inside the magnetic layers. However, since the study of the CPP- 
it is well known that, in this geometry, the spin polarization of the current is 
strongly affected by the spin accumulation which is controlled by spin relaxation effects in 
the whole structure. To correctly describe the electronic transport, diffusive models that 
also take into account this relaxation must be considered^. Indeed recent spin transfer 
models do take into account these diffusive transport aspects^ 8 * 12 * 13 * 14 * 1 -^. 

The first experimental evidence of diffusive effects on the STT was given by Urazhdin 
et alW^ and Alhajdarwish et alW^. They demonstrated the influence of the spin dependent 
scattering by magnetic impurities on the STT features in the CIMS regime at low magnetic 
field : by doping the polarizing layer with impurities, the amplitude and even the sign of 
the torque was changed. More recently, we demonstrated experimentally the influence of 
the spin relaxation effect on the STT^. We showed that by acting on the distribution of 
spin relaxation in the structure through different spin diffusion lengths in the two mag- 
netic layers, the STT and thus the dynamics of magnetization induced by this torque can 
be strongly modified compared to standard structures. This was made by studying the 
magnetization dynamics induced by the STT in a Co(8 nm)/Cu/Py(8 nm) (Py=NigoFe2o) 
nanopillar through both static and high frequency magneto-transport measurements. In 
such a structure, an unusual dependence of the STT with the angle <p between both magne- 
tization is predicted in the frame of diffusive models of spin transfer ^ 3 * 19 * 20 * (see Fig. [T]). This 
so called Wavy Angular Dependence of the Spin Transfer Torque (WAD-STT) arises from 
the modification of the spin accumulation profile in the structure caused by the different 
ratios thickness /spin diffusion length of the magnetic layers. It is characterized by a change 
of sign of the STT between (parallel configuration (P)) and it (antiparallel configuration 
(AP)). This angular dependence modifies the stability of the static states P and AP: for 
one polarity of the current (negative current in our convention), both states are stabilized 
by the STT whereas they are destabilized for the opposite polarity (positive current). In the 
latter case, it is predicted that the magnetization precesses even at zero external magnetic 
field. These predictions were validated experimentally by the present authors by measuring 
the microwave emission caused by the steady precession of magnetization induced by the 
STT at zero (or low) applied fielcP^. 

In this manuscript, we present a complete characterization and analysis of the magneti- 
zation dynamics in the presence of a WAD-STT. The results are compared to the prediction 
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of the magnetization dynamics obtained by single domain magnetic simulations (macrospin 
approximation). We show that the main static and dynamical experimental features such 
as current-field phase diagram, variation of resistance with current, variation of frequency 
with current and field can be reproduced by simulation, at least qualitatively and for some 
features quantitatively. Therefore, the present work validates the diffusive models of the 



The paper is organized as follows. We first described the predictions of the macrospin 
simulations taking into account the WAD-STT. In a second section, we present the results 
of low field static and high frequency measurements obtained on a Co(8 nm)/Cu/Py(8 nm) 
nanopillar in which a WAD-STT is predicted. These results are compared with the predic- 
tions of the macrospin simulations and discussed in a third part . 

2. MACROSPIN SIMULATIONS WITH A WAVY ANGULAR DEPEND ANCE 
OF THE SPIN TRANSFER TORQUE (WAD-STT) 

2.1. Method 

In this section, we present the results of the simulations for the magnetization dynamics 
of the Py free layer of a Co(8 nm)/Cu/Py(8 nm) trilayer, on the basis of the predicted 
wavy angular dependence of the torque (labelled as WAD-STT). These calculations have 
been carried out in the macrospin approximation, i. e by assuming the magnetization to be 
uniform during its motion. The magnetization of the Py free layer is described by a unit 
vector m = m x u x + m y u y + m z u z and makes an angle ip with the magnetization of the 
Co layer which is supposed fixed and aligned along the x axis. As in experiments, the Py 
free layer has got an elliptical shape of dimensions 100x155 nm 2 and a thickness tp y of 8 
nm. This results in an in-plane uniaxial anisotropy field H an = H an m x u x aligned along the 
direction of the long axis u x of the ellipse and a demagnetizing field Hd = —Hdin z u z , with u z 
perpendicular to the plane of the layer. The external field H app = H app u x is applied in the 
plane of the layer and aligned along the direction u x of the uniaxial anisotropy field. The 
time dependent trajectory rh(t) is found by solving the modified Landau-Lifschitz-Gilbert 
(LLG) equation including the STT r: 



spin transfer^™!! 



dm 
~dt 



—70 rh x Hf + a rh x 



dm 70 



dt HqM s V 



T 



(1) 
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Here 70 is the absolute value of the gyromagnetic ratio, a the parameter of the Gilbert 
damping, M s the saturation magnetization, V the volume of the Py nanomagnet and Hf 
an effective field which includes H app , H an and H^. In order to characterize the influence of 
thermal activation on the magnetization dynamics, the simulations were carried out both at 
T=0 K and T=300 K. Thermal effects are simulated by introducing a randomly fluctuating 
field w ith hqHt = ^/2k B Ta/VM s -fAtH a i, where H at is a random gaussian field 

with < H a i >= and < >= 1, kg is the Boltzmann's constant and At the integration 
step. The magnetic parameters used for the simulations are described in appendix |A| with 
computational details. The spin transfer term is defined as r = P(y)p m x m x u x , I is 
the current (defined as positive when the electrons flow from the fixed Co layer to the Py 
free layer), e is the absolute value of the electron charge. The polarization factor P(<p) is 
derived from the angular dependence of the torque calculated from the Barnas-Fert model 
(see Fig. [TJ. 

The displayed phase diagrams are constructed by calculating the mean values over the 
integration time of a mx = ■J < in* — < vn x > 2 > after relaxation (a mx gives an evaluation of 
the amplitude of precession). These two quantities allow to determine the borders between 
the different static and dynamical states. All presented simulations have been carried out for 
a positive current. No excitations are predicted by the simulations for a negative current!^. 

2.2. Results 

2.2.1. T=0 K 

In Fig. |2^a) and (b), we show respectively the variation of the frequency and the resistance 
as a function of the current for zero applied field (black curves). Starting from an initial 
P configuration and increasing current (plain line), the magnetization goes into sustained 
precession around the in plane effective field (in plane precession mode IPp) for a threshold 
current l\ (black point at 1=4.8 mA) with a precession frequency (f=2.28 GHz) close to the 
Kittel's frequency of the FMR small angles precession (f=2.23 GHz). When the current is 
increased above l\, the angle of precession increases but the frequency decreases ("red shift 
regime"). This behavior, i.e. the decrease of frequency with increasing current is a nonlinear 
effect due to the dependence of the frequency on the precession amplitude^ 1 ™ 1 . In this 
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mode, only a small increase of the resistance with increasing current is predicted, although 
a large angle precession can be reached. For a threshold current J? + (orange point at 
1=12 mA), the magnetization dynamics changes abruptly to an out-of-plane precession mode 
(labelled as OPp ). This transition is associated with a jump in frequency and a large increase 
of the resistance. In this OPp regime, the magnetization precesses around the demagnetizing 
field and the frequency is set by the mean out-of-plane component of magnetization < m z >. 
For increasing current, the trajectory goes away from the plane of the layer (higher < m z >) 
and the frequency increases ("blue shift regime"). For higher current (I> 13 mA, not shown), 
the magnetization precesses faster (up to 20 GHz) around circular trajectories of decreasing 
radius while current is increased. These trajectories converge for even higher current to a 
quasi-static out-of-plane state OPs- As shown in Fig. |2](a-b), the transition between the 
IPp and the OPp precession modes is irreversible): the critical current for this transition is 
higher for increasing current (/,? + , orange point and plain line) than for decreasing current 
{Ic~ , blue point and dashed line), which results in hysteretic R(I) and f(I) curves. 

In addition, we display in Fig. [2]^a) the variation of the frequency as a function of the 
current for a higher negative applied field H = -20 Oe (red curve). In the IPp mode, the 
black curve (H=0 Oe) is above the red curve (H=-20 Oe): the frequency increases with 
field. On the contrary, in the OPp mode, the black curve is below the red curve and the 
frequency decreases with field: the in-plane field pulls magnetization in the plane of the 
layer, decreasing < m z > and therefore the frequency. 

We present in Fig. [2|c) a calculated current-field dynamical phase diagram for increasing 
current starting from an initial P configuration. The amplitude of precession (a mx ) is plotted 
in a color scale. Several regimes can be defined. First, at high positive field (H > +25 Oe) 
(see the orange dashed line in Fig. [2?c)), only the IPp mode is observed for I\ < I < 1%. 
For I > 1%, the magnetization goes into the quasi-static out-of-plane precession state OPs- 
In a second regime at low field i.e. for fields smaller than the coercive field (see the white 
dotted line on Fig.|2](c)), the out-of-plane precession mode OPp appears for I above Jf + with 
I\ < Il + < I < 1%. The amplitude of precession (red color) and therefore the emitted power 
is maximum in this mode. The aforementionned hysteretic behavior appears also on the 
phase diagram: the OPp /IPp frontier is not the same for increasing current (colored frontier, 
7? + ) and decreasing current (dashed lines, /?~). At higher current and field (H > +25 Oe), 
the transition to the OPs state is also irreversible: for increasing current (colored frontiers), 
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magnetization switches from IPp to OPs and for decreasing current from OPs to OPp 
(dotted lines). 

Finally, for even a lower applied field (H< —50 Oe), before reaching the OPp mode, 
the magnetization first precesses in-plane around the AP state (mode denoted as IPap). 
Note that in this mode, the frequency increases with current and decreases with (positive) 
field. Starting from an initial AP state (H< —90 Oe), this precession state IPap is the first 
observed when current is increased from zero. 

To summarize, the magnetization dynamics in a WAD-STT structure is very different 
from what is predicted for a standard angular dependence structure: 

• For H > H c , where H c is the coercive field of the free layer, the precession modes are 
predicted for a positive current for a WAD-STT structure in our current convention, 
whereas they are obtained for a negative current for a standard structure. 

• For H < H c , in the standard case, only small angle IPp precessions around the in-plane 
effective field are predicted in a very narrow range of current before the magnetization 
switches to the other stable static state (P or AP) due to the increase of the precession 
amplitude with currentP^Hl. In the case of a WAD-STT, both static states (P or AP) 
are unstable and this in-plane-precession mode is present on a much larger range of 
current. In a similar manner to what is observed for a standard angular dependence 
at high field, large angle in-plane trajectories bifurcate into out-of-plane trajectories 
at higher currentpSES. One can note that this large angle out-of-plane precession is 
predicted even in the absence of any external field. It is therefore possible in this 
structure to emit microwaves at zero field, only by injecting a current in the structure. 

2.2.2. T=300 K 

The main effect of temperature is the almost complete disappearance of the irreversibility 
associated with the transition from the IPp to the OPp precession mode. The thermal energy 
allows the magnetization to go across the energy barrier that separates these two states. 
This appears on the phase diagram of Fig. [3]^c): the frontiers between the two states are 
quasi-identical for increasing (colored frontier) or decreasing current (dashed lines), except 
at large applied field (H > 60 Oe). As a consequence, this reversible behavior leads to an 
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expansion of the current-field region in which the OPp precession mode occurs: for a given 
value of the applied field, the critical current for the transition from IPp to OPp is lower 
at T=300 K^. The reversibility of this transition results in non-hysteretic f(I) and R(I) 
curves as shown in Fig. j3^a-b))^l The effect of thermal energy appears in particular in 
the transition zone between the two states and especially for currents close to the minimum 
in frequency (I~ 9.1 mA). In this zone, the thermal energy makes the system unstable 
and telegraph noise with a characteristic time of about 10 ns is observed between the state 
IPp and the two degenerated out-of-plane trajectories which are symetrical around the (xy) 
plane. To conclude this section, the temperature does not modify the main characteristics 
of the dynamics in a WAD-STT structure, such as a steady precession at zero applied field. 

3. EXPERIMENTAL RESULTS AND COMPARISON WITH MACROSPIN SIM- 
ULATIONS 

We present a detailed experimental characterization of the magnetization dynamics in- 
duced by STT on a Co(8 nm)/Cu(10 nm)/Py(8 nm) nanopillar in which a WAD-STT 
is predicted. This dynamics was studied by measurements of the high frequency voltage 
oscillations generated by the magnetization precession with a spectrum analyzer after am- 
plification. To directly compare the experimental results with simulations, we consider for 
experiments an effective external field H e ff that takes into account the contribution from 
the dipolar field exerted on the Py free layer (Hdi P ~ —43 Oe). 

3.1. Initial configuration : P 

In Fig. |4]^a), we display the microwave spectra measured for several values of injected 
current and a fixed effective applied field H e ff = — 62 Oe . The corresponding variation of 
the frequency with current at this field is plotted in Fig. [4](b) . The analysis of the frequency 
dependence with current allows the definition of three different dynamical regimes. At low 
current (4 mA^ 7^9 mA), we observe that the frequency increases with current i.e. 
a blue shift with current (regime 1). Then it stays approximately constant for 9 mA^ 
/ ^ 10 mA) (regime 2). At higher current (/ > 10 mA), a slight decrease of frequency 
with current is observed, i.e. a red shift with current (regime 3)^. Three corresponding 
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behaviors are observed in the frequency dependence with field as illustrated in Fig. |3|c) 
by the microwave spectra measured for a constant current of 10 mA and several values of 
applied field. The corresponding frequency are plotted Fig. md)). In the low field range (- 
47 Oe^ H e ff ^ 2 Oe), the frequency decreases with field (regime 1). At higher negative field 
(-62 Oe^ H e ff ^ —47 Oe), the frequency remains approximately constant (f~ 3.5 GHz) with 
field (regime 2). Then in a narrow experimental window for H e ff very close to the switching 
field, the frequency decreases with field (regime 3). 

In Fig.^e), we present a current-field dynamical phase with the microwave power plotted 
in color scale. As predicted by the model, the microwave power is emitted for a positive 
current and at low field {H e ff < 30 Oe), i.e smaller than the coercive field, and in particular 
at zero effective field. This behaviour is very different from the one typically observed in 
standard structures, where steady precession associated with microwave emission is observed 
at high field, i.e higher than the coercive field, and for a negative currentP^. In addition, 
we plot on the diagram in plain lines the frontiers between the different dynamical regimes 
presented above^l In regime 1, the frequency increases with current and decreases with 
field. It is present in the largest part of the region of the phase diagram where a power 
emission is observed. This behaviour is indeed the only one observed at low current (from 
1^ 4 mA) and close to zero effective field. The regime 2 appears at low applied field and 
for higher currents. It is characterized by a frequency that remains approximately constant 
with the current and also with the applied field. Finally, the regime 3 in which the frequency 
slightly decreases with current (and increases with field) is observed in the top-left corner 
of the diagram corresponding to large current and applied fields close to the switching field 
(H eff « -90 Oe). 

We can directly compare the experimental results to the simulated dynamical phase dia- 
gram shown in Fig. [4](f ) . The current-field region in which an out-of-plane precession OPp 
mode is predicted (red color in Fig. |4|e)) coincides with the zone in which the high fre- 
quency power is measured (dynamical regime 1, 2 and 3. Note that the field dependence 
of the experimental critical currents is in excellent agreement with the one associated with 
the IPp/OPp transition in the simulations. In addition no microwave power was measured 
in region where an IPp mode is predicted . We believe that, most probably, the weak 
microwave emission associated with these IPp mode in a WAD-STT structure is too small 
to be measured with our experimental setup (this point will be discussed in the section El). 
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Finally, the predicted IPap precession mode in the lower left corner of the diagram char- 
acterized by a red shift regime with current and moderate precession amplitude was not 
observed experimentally. 

3.2. Initial configuration : AP 

We can go deeper into the compared analysis of experiments and simulations by looking at 
some additional features observed when one starts from the antiparallel AP configuration. In 
Fig. [5(a) and (b), we show respectively the experimental and the simulated phase diagrams 
for increasing current starting from an initial AP configuration. The experimental diagram 
was constructed by measuring R(I) curves for different values of applied field and increasing 
the current starting from the high resistance AP configuration. The resistance level is 
plotted in colour scale in Fig. [5(a) and the corresponding R(I) curves are plotted Fig. [5(c). 
Experimentally, starting from the AP configuration and above a positive threshold current, 
the magnetization switches sharply from the AP to P configuration. Then for a larger 
current, an increase of the resistance is observed due to the onset of a steady precession 
mode. This behavior is well reproduced by the macrospin simulation (see Fig. ||(b)). In 
particular, the critical currents associated with the magnetization switching from AP to P 
state are well reproduced, as well as their dependence with the applied magnetic field. 

Now we compare the experimental (Fig. [5(c)) and simulated R(I) (Fig. [5(d)) curves ob- 
tained for increasing current starting from an initial AP configuration. A good agreement is 
obtained (see for example the critical currents or the variation of resistance). In particular, 
the simulated increase of resistance associated with the transition from the IPp to the OPp 
state is reversible as is the case in experiments. If one calculates the differential resistance 
dV/dI=R+IdR/dI using the simulated R(I) curves, one observes that this transition is as- 
sociated with a peak in the dV/dI(I) curves. This reproduces the experimental observations 
where such peaks, that characterize the reversibility of the transition, have been always ob- 
served at the onset of microwave emission (see Ref.n®). In addition, for both simulations and 
experiments, no significant variation of the resistance is observed for the in-plane preces- 
sion IPp mode. However, the simulation overestimates the stability of the IPap precession 
mode around the AP configuration at the expense of the P state (see IPap zone in Fig. plb) 




and Fig. §d)). This can be explained by an underestimation of the effect of temperature 
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caused by a sweeping ramp 10 5 faster in simulations than in experiments. With a slower 
current ramp, the I Pap mode should be more unstable and thus the reversal toward the P 
state should be favored. This may also explain why we did not measure microwave power 
associated with this dynamical precession mode. 

3.3. Frequency of the precession 

We compare in Fig. |6](a) the experimental and calculated f(I) curves for several values 
of applied field. A good agreement is obtained for fields close to the zero effective field 
(H e // = 2 and 10 Oe) corresponding to the regime 1 in the phase diagram. In particular, 
the blue shift with current and the decrease of frequency with field observed experimentally 
is well reproduced by the simulation. Simulations and experiments significantly differ for 
higher negative field (H e // = —45 Oe). In particular, the saturation regime (regime 2) 
characterized by a frequency approximately constant with current is not reproduced by 
the simulation. However, the experimental slope of frequency with current df/dl before 
saturation (H e// = -45 Oe (7 mA^ I ^ 9 mA), H e// = 2 Oe and H e// = 10 Oe) is well 
reproduced by the simulations. This seems to indicate that the regime before saturation 
(regime 1) is closer to an homogeneous precession mode than the regimes observed for 
higher current (regime 2 and 3). 

The disagreement simulation/experiments of the f(I) curve at high negative field reflects 
a decrease of frequency with field lower in simulations than in experiments. This appears 
on Fig. |6](b-c), where the experimental and simulated dependence of the frequency with 
magnetic field are compared for 1= 10 and 11 mA. However, the general shape of the exper- 
imental f(H) curve in the red shift regime for higher field values is qualitatively reproduced 
by the simulation. 

4. DISCUSSION 

The macrospin simulations well reproduce in overall the main experimental features: the 
region in current-field diagram with larger power emission and increased resistance, the 
frontiers between the main different static and dynamical states, and the dc resistance of 
these different states. This agreement is quite satisfactory if we consider the crude macrospin 
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approximation that was used. 

The direct comparison of the experimental and simulated dynamical phase diagrams 
clearly suggests that the observed excitations are associated with an out-of-plane precession 
mode. This is confirmed by the blue shift in current and the red shift in field observed in the 
dynamical regime 1, signature of this dynamical mode. However, the macrospin simulation 
do not reproduce the saturation regimes observed at low applied field and for a large current 
(regime 2 and 3) (see Fig. [6^ a)). In this range of parameters, we believe that these regimes are 
probably associated to an inhomogeneous distribution of magnetization. It is worth noting 
that the only two other groups who observed a blue shift regime for an in-plane applied field, 
did also observe this saturation (or decrease) of the frequency for larger applied current (see 
Fig. 3, RefP^ and RefP^), suggesting the intrinsic character of this dynamical feature. On 
the basis of micromagnetic simulations®', it has been recently suggested that this behavior 
can be attributed to the onset of inhomogeneous dynamical modes characterized by a mean 
out-of-plane component of magnetization < m z > lower than for the uniform mode. The 
onset of such non standard modes leads to a decrease of the large magnetostatic energy 
generated by the strong demagnetizing field when the magnetization goes out of the film 
plane. The frequency of precession being proportional to < m z >, the decrease of < m z > 
results in a decrease of the frequency as compared to the macrospin prediction and thus 
explains this saturation phenomenon. The in-plane Oersted field may also play a role: in 
the case of an out-of-plane precession mode, it tends to pull the magnetization towards the 
film plane. This also leads to a decrease of the precession frequency^. 

Another point to discuss is, as already mentioned, the non-observation of a red shift 
regime associated with the IPp mode in our experiments. This regime was not ob- 
served experimentally in the Co(8 nm)/Cu/Py(8 nm) nanopillar described in this pa- 
per. However, we could observe it in a nanopillar having a different WAD-STT structure 
Co(4 nm)/Cu/Py(8 nm). In this nanopillar, before the onset of the higher power blue shift- 
ing microwave peaks, signals of very low amplitude presenting a clear red shift behaviour 
have been observed. In standard nanopillars, this IPp mode corresponding to a large pre- 
cession angle goes with a large emitted power^S This difference between standard and 
WAD-STT samples can be explained by a weaker slope of the angular dependence of the 
GMR for angles around cp = in the case of a WAD-STT 33 . This results in a weak variation 
of the resistance (and therefore a very weak power) for in-plane oscillations of magnetization 
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even for large angles. In our calculations, we accounted for this specific angular dependence 
of the GMR by using a normalized resistance r(cp) = sin 2 (<^/2)/(l + x cos 2 (<^/2)) with an 
asymmetry factor^ ^ = 7.7. As can be seen on Fig. [3^b) or reffig4(d) , this results in a low 
increase of resistance of the device in the IPp mode. To evaluate the influence of the angular 
dependence of the GMR on the emitted power, we calculated the maximum emitted power 



for this IPp mode using Eq. Al with x — 7.7, x = 2 (measured experimentally by Urazhdin 
et al. in a standard structure Py(6 nm)/Cu/Py(12 nm), and x = 0) corresponding to the 
simple angular dependence AR = (1 — cos6 l )/2. The simulated power in the precession mode 
for x — 7.7 is about 10 times lower than for x — 2 and 25 times lower than for x = 0. If 
this assumption is valid, then it allows to explain that in our experimental configuration, 
the measured signals in this IPp precession mode were too weak to be detected. 

Finally, our simulations allow the calculation of the output power in the OPp mode: it 
is about 3.9 pW/mA 2 , should compare to the maximum measured output power of about 
4.5xl0~2 pW/mA 2 , i.e approximately a factor 80 smaller. This discrepancy may arise 
from several factors. First, due to the impedance mismatch and the attenuation of the 
signal in the high frequency line and the variation of the amplification gain with frequency, 
the measured signal (after taking into account the amplification) is only a fraction of the 
actual emitted signal. Second, the fact that the OPp mode predicted by the macrospin 
simulation has been rarely observed experimentally^ indicates the macrospin approximation 
does not describe properly the dynamics in this precession mode. Recently, micromagnetic 
calculations by Berkov et a/J 34 * 35 ^ have shown that this dynamical OPp mode is sensitive to 
the magnetization homogeneity. Factors favoring the inhomogeneity such as temperature or 
the Oersted field decrease the output power in this precession mode. Third, the possible 
presence of ferromagnetic oxydes (NiO) on the sidewalls of the pillar due to air exposure 
during the fabrication process®' or the stronger Oersted field on the edges may force the 
magnetization to lie in the plane on the edges. In this case, only a part of the magnetization 
in the center might therefore precess in the out-of-plane direction, the resulting output power 
being weak. 
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5. CONCLUSION 



As a conclusion, we have simulated the magnetization dynamics in the macrospin ap- 
proximation taking into account the WAD-STT predicted for a Co(8 nm)/Cu/Py(8 nm) 
nanopillar by the Barnas-Fert modeP^^. At low field, the magnetization dynamics is 
strongly modified as compared to a standard angular dependence of the torque. When the 
current is increased, at low and even zero applied field, one observes successively an in- 
plane precession mode (IPp) in which frequency decreases with current ("red shift regime") 
and an out-of-plane precession mode (OPp), in which the frequency increases with current 
("blue shift regime"). The transition between both regimes is associated with an increase 
of the dc resistance and a strong increase of the output power. We have presented some 
experimental results that confirm this zero and low field steady precession. In the main 
part of the current-field phase diagram in the low field range, the frequency increases with 
current (blue shift regime) and decreases with field. In the remaining part of the diagram 
(large negative applied field and high applied current), a regime in which the frequency 
is approximately constant (or decreases slightly) with current is observed. The macrospin 
simulations reproduce the main observed experimental features (phase diagram, R(I) curves, 
dependence of frequency with current and field) and suggest that the dynamical excitations 
observed experimentally are associated with an out-of-plane precession mode. However the 
saturation regimes observed at lower field and higher current are not reproduced by the 
simulation and may reveal the onset of inhomogeneous dynamical excitations. Following 
the first experimental confirmation^ of the STT induced zero field steady precession caused 
by the WAD-STT, this detailed comparison of the static and dynamical results between 
models and experiments and the observed general good agreement give an additional con- 
firmation of the predicted WAD-STT by diffusive models of spin transfer. This possibility 
of engineering the angular dependence of the torque with different spin diffusion lengths in 
the magnetic layers underlines the role of the spin accumulation and the spin relaxation 
effects on the STT. By playing on the distribution of spin relaxation in the structure, it is 
therefore possible to strongly modify the STT and the magnetization dynamics induced by 
the torque. 
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APPENDIX A: PARAMETERS USED FOR THE MACROSPIN SIMULATION 

For the macrospin simulation a gyromagnetic ratio 70 = 2.2. 10 5 m/(A.s) has been used. 
The saturation magnetization fioMs = 0, 87 T has been deduced from ferromagnetic reso- 
nance experiments carried out on Cu(6nm)/Py(7nm)/Cu(6nm) thin films. The anisotropy 
fields Ha and H an can be expressed as a function of demagnetizing factor N x , N y et N z as- 
sociated with the shape of the nanomagnet with Hd ~ M S (N Z — N y ) et H an « M s (N y — N x ). 
If we approximate the shape of the nanomagnet with an ellipsoid with axes 150 nm, 105 nm 
and 8 nm, close to the lateral dimensions measured by SEM within measurement uncertainty, 
one obtainsS3Nx=0.047, Ny=0.063, Nz=0.89 leading to fi H d « 0.7 T and (i H an « 0.014 T. 
This value of H an is close to the anisotropy field HoH' an « 0.0145 T one can estimate from 
the room temperature experimental coercive field H c m . For a direct comparaison between 
experimental and theoretical phase diagrams, the uniaxial anisotropy field H an was adjusted 
to obtain a simulated coercive field equal to the experimental one (hqH c = 0.009 T). For 
simulation at T=0 K, we used therefore ^oH an = 0.009 T. At T=300 K and with the pa- 
rameters and waiting time used in the simulation, one obtains ^H an = 0.012 T. A Gilbert 
damping parameter a = 0.018 has been used. Finally, currents are deduced from current 
densities using a lateral surface of the nanopillar of A=1.38xl0~ 14 m 2 , measured experi- 
mentally by SEM. Since the angular dependence of the GMR in a WAD-STT structure is 
predicted 33 * 38 ^ to deviate significantly from the commonly used dependence sin 2 (<^/2), the 
normalized magnetoresistance r = (R — Rp)J (Rap — Rp) was calculated using the following 
relation™! 

r(<p) = sin 2 (^/2)/(l + xcos 2 (v?/2)) ( A l) 
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We have used x = 7.7 that was derived experimentally by Urazhdin et a/P^ in a non 
symmetrical pillar composed of Py(6 nm)/Cu(10 nm)/Py(1.5 nm) for which a WAD-STT 
is also predicted. 

The precession frequency was derived from the higher amplitude lowest order peak in the 
Fourier transform spectrum of the m y component. In the case, very common experimentally, 
of a slight misalignment between the magnetization of the polarizing layer and the easy axis 
of the free layer, this frequency corresponds actually to the higher amplitude lowest order 
frequency measured experimentally^^!!. The output power P corresponding to the simulated 
trajectories is deduced from the reduced resistance r using P/I 2 = AR 2 . /Z c . < (r— < r > 
) 2 >, with AR exp = Rap — Rp = 51 mfl the experimental static variation of resistance due 
to the GMR, and Z c = 50 Q the characteristic impedance of the high frequency line. < ... > 
indicates the mean value of 40 ns after relaxation of the magnetization. 

The LLG equation (jlj has been solved using a fourth order Runger-Kutta algorithm 
with a time step of 1 ps. For simulation at T=0 K, to let magnetization reach a stationnary 
state, a 100 ns relaxation time has been used, then trajectories are saved over a 40 ns. 
These trajectories are used for the fast Fourier transform calculations and the deduction 
of the magnetoresistance. For simulation at T=300 K, after application of the current, a 
relaxation time of 30 ns has been used and simulations have been carried out by sweeping 
current with a 10 5 mA/s ramp. The plotted frequency and resistance vs current curves are 
averaged over 5 realizations or more. 
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The coercive field H c at a temperature T can be expressed as a function of the uni- 
axial anisotropy field H' a ^ H c (H' an ) = H an (l - [k B T/E k ln(/ r m / In 2)] 1 / 2 ) with E k = 
l/2(j,oM s VH' an , /o the attempt frequency of the order of 1 GHz and r m the time measure- 
ment of the order of Is. By solving this equation, using experimental parameters and assuming 
/ =1 GHz, one finds H' an m 145 Oe. 

One can note however that the frequency derived from the m x and the m y component are the 
same for an out-of-plane precession mode. 
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1/2 

Fig1, 
Boulle et al. 

FIG. 1: Magnitude of the spin transfer torque r on the free Py layer of a Au(infinite) /Cu(5 
nm)/Py(8 nm)/Cu(10 nm)/Co(8 nm)/Ta(10 nm)/Cu(infinite) multilayer as a function of the angle 
(p between the magnetizations of the free Py and fixed Co layers, calculated with the Barnas-Fert 
model 13 (r(ip) = — P(ip) sin iplh/(2e)). Electrons flow from the Co layer to the Py layer for a 
positive current I. 



21 




4 6 8 10 12 
Current (mA) 

Fig. 2, Boulle etai. 



FIG. 2: T=0 K. (a) Frequency of precession as a function of current for several applied field for 
increasing (plain line) and decreasing (dashed line) current, (b) Normalized resistance as function 
of current for H=0 for increasing (plain line) and decreasing (dashed line) current, (c) Dynamical 
phase diagram for increasing current (o mx in color scale) and decreasing current (dashed and dotted 
frontiers). For decreasing current, the dotted white line corresponds to the frontier of the OPs 
with the OPp state and the dashed black line to the frontier of the OPp state with the IPp state. 
Plain black lines define the frontiers of the P and AP state. Static and dynamical states are defined 
in the text. 
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FIG. 3: T=300 K. (a-b) Frequency of precession (a) and normalized resistance (b) as a function 
of current for increasing (upward triangles) and decreasing (downward triangles) current for H=0. 
(c) Dynamical phase diagram for increasing current {a mx in color scale) and decreasing current 
(dashed and dotted frontiers). 
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FIG. 4: (a): Microwave power spectra for several values of injected current and H e ff = — 63 Oe. 
(b) Frequency associated to the peaks in the microwave spectra as a function of current for H e ff = 
—63 Oe. (c): Microwave power spectra for several values of H e ff and 1=10 mA. (d) Frequency 
associated to the peaks in the microwave spectra as a function of the H e ff for 1=10 mA. H e ff is 
defined as ~H. e ff = H app + Hfc p with H^ p = —43 Oe. In (b) and (d), error bars correspond to the line 
width of the peak, (e) Experimental dynamical current-field phase diagram for increasing current 
starting from an initial P state. The microwave power is plotted in color scale. The dynamical 
regimes denoted as 1, 2 and 3 are defined in the text, (f) Simulated current- field phase diagram for 
increasing current starting from an initial P state. The amplitude of oscillation (<7 m;r ) is plotted 
in color scale. 
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FIG. 5: (a) Experimental phase diagram for increasing current starting from an initial AP state. 
The normalized resistance of the sample AR is plotted in color scale, (b) Corresponding simulated 
phase diagram (normalized resistance in color scale), (c): Measured resistance as a function of 
current for several value of H e ff starting from an initial AP state for successively increasing and 
decreasing current, (d): Corresponding simulated normalized resistance vs current curves. The 
experimental normalized resistance AR is obtained by substracting from the experimental R versus 
/ curves a reference curve to remove the changes in resistance due to Joule heating. 



25 



N 



3- 



>* 
o 
c 

O Q 



(a) 

10 Oe I 

2 0e | I 

" -45 Oe ■ 

■ 

■ / 


l i ¥ m M 

/ w/ * 

// 1 


i.i. 


■ Experiment 

— Simulation 
i . i 



8 10 
Current (mA) 



12 




(C) 



— ■—11 mA 
—■-10 mA 



4. 

3. 

3. 
2, 
2. 



N 



c 

C) 

cr 



-40 -20 

H e„<Oe) 



-80 -60 -40 -20 

H e „( 0e > 



20 



Experiment Simulation 
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FIG. 6: (a) Simulated (plain lines) and experimental (square points) variation of frequency with 
current for several values of H e ff. (b-c): Simulated (b) and experimental (c) variation of frequency 
with H e ff for 1=10 mA and 1= 11 mA. 
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